// project home: cp3-cg.googlecode.com

#ifndef __MVP_H__
#define __MVP_H__
#include <Rtypes.h>
#include <vector>
using namespace std;

vector<Double_t> matrixvectorproduct(vector<Double_t> *A, vector<Int_t> *JA, vector<Int_t> *IA, vector<Double_t> *x,bool debug=false);

Double_t skalarp(vector<Double_t> *x, vector<Double_t> *y);


// the returned vector is x-y
vector<Double_t> vectorsub(vector<Double_t> *x, vector<Double_t> *y);

// the returned vector is x+y
vector<Double_t> vectoradd(vector<Double_t> *x, vector<Double_t> *y);


// the returned vector alpha*x
vector<Double_t> VectorMultiS(vector<Double_t> *x, Double_t alpha);




#endif